Vortex core deformation and stepper motor behavior in a superconducting ratchet 
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We investigated experimentally the frequency dependence of a superconducting vortex ratchet 
effect by means of electrical transport measurements and modeled it theoretically using the time 
dependent Ginzburg-Landau formalism. We demonstrate that the high frequency vortex behavior 
can be described as a discrete motion of a particle in a periodic potential, i.e. the so called stepper 
motor behavior. Strikingly, in the more conventional low frequency response a transition takes place 
from an Abrikosov vortex rectifier to a phase slip line rectifier. This transition is characterized by 
a strong increase in the rectified voltage and the appearance of a pronounced hysteretic behavior. 



Since the early 90's, inspired by the pioneering work of 
Magnasco [1 , a strong theoretical effort has been devoted 
to modeling complex biological ratchet systems using the 
one-dimensional overdamped Langevin equation of mo- 
tion for a single particle |2l [Sj. In this approximation, 
the possibility of a directional motion was demonstrated 
for a Brownian particle, affected by an external time- 
correlated driving force and a spatially periodic potential, 
which lacks reflection symmetry. A common prediction 
of these models is discretization of the particle dynamics 
under a periodic driving force: the average velocity of 
the particles can only acquire discrete values correspond- 
ing to an integer number of unit cells traveled by the 
particles during each cycle of the drive. Soon after these 
theoretical predictions were published, the main concepts 
have been successfully applied in experimental research 
both at microscopic and macroscopic scales |34 i4j . A par- 
ticularly ideal playground with a high level of flexibility 
allowing to explore the influence of particle density, parti- 
cle size, and to cover the whole space of parameters of the 
ratchet potentials is provided by vortices in a type-II su- 
perconducting film [5]. Here, the particle is a quantized 
flux bundle (vortex), which is induced by an external 
magnetic field and which has a temperature-dependent 
size. These vortices can be driven by an applied electrical 
current flowing through the superconductor whereas the 
periodic asymmetric potential to trap the vortices can be 
generated by using modern lithographic techniques [6l [7] . 

However, two important features inherent to vortex 
ratchets have not been properly taken into account by the 
current theoretical models. Firstly, the repulsive interac- 
tion between vortices makes their dynamics more com- 
plex. Indeed, the competition between the asymmetric 
pinning and the vortex- vortex repulsion leads to reversals 
of the preferential flow direction of vortices as their den- 
sity increases [H [9] . Secondly, the one-to-one mapping of 
a vortex to a particle breaks down at high vortex veloci- 
ties where vortex displacements take place at shorter time 
scales than the healing time to recover the superconduct- 
ing order parameter along the vortex trajectory. As the 



vortex velocity is progressively increased, this latter ef- 
fect leads first to a dynamic reconfiguration of the vortex 
lattice, followed by the elongated- vortex motion which 
is much faster than that of the Abrikosov vortices. The 
driven elongated vortices are eventually transformed into 
phase slip lines [lOl [11] . This bring us to the formulation 
of the main unsolved problems that will be addressed in 
the present work, (i) how do the internal degrees of free- 
dom of a vortex influence the rectification properties of 
vortex ratchets? (ii) is it possible to realize the theoret- 
ically predicted but so far experimentally not observed 
vortex stepper motor behavior? 

Experimental details and theoretical model. -In order to 
address these important issues we have investigated a su- 
perconducting vortex rectifier consisting of a 33 nm thick 
Al bridge patterned with a periodic array of asymmet- 
ric pinning sites with a unit cell consisting of two square 
antidots of different sizes (1200 nm and 600 nm), sepa- 
rated by 400 nm distance [See Fig. [ija)]. The supercon- 
ducting critical temperature is Tc= 1.445 K. The coher- 
ence length and the penetration depth are ^(0) ~ 73 nm 
and A(0) ~ 261 nm, respectively. When applying an ac 
current along the bridge (x-axis), the vortices, induced 
by a magnetic field perpendicular to the plane of the 
bridge (2;-axis), are subjected to a driving force in the 
^/-direction. The resulting motion of vortices leads to dis- 
sipation, which is probed by measuring the voltage drop 
along the bridge. The superconducting vortex rectifier 
has been investigated both experimentally in a ^He cryo- 
stat and theoretically using time-dependent Ginzburg- 
Landau (TDGL) simulations. The previously used the- 
oretical approaches (Langevin dynamics, Fokker-Planck 
equations, molecular dynamics), which are not specific 
for superconducting systems, deal with motion of rigid 
entities embedded in a medium, while vortices in a super- 
conductor correspond to "soft" superconducting current 
loops which can be deformed by applying exteranl drive 
suchas the Lorentz force. Unlike these approaches, the 
TDGL formalism [121 [13] treats the dynamics of the con- 
densate as a whole, in terms of the time-dependent com- 
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plex order parameter ^/^(r, t). This treatment contains all 
necessary ingredients, including vortex deformability, to 
reproduce and explain the experimental data. 
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FIG. 1. (a) Atomic Force Microscopy image of a 5 x 5 
/xm area of the sample, (b) 1-01^ and (c) Re calculated 
at if = 2Hi using the TDGL model, with Hi the magnetic 
field generating exactly one single-quantum vortex in each 
unit cell. Experimentally obtained: (d) resistance versus field 
measured at 1.2 K, and with dc current of 100 /iA. (e) The 
dc voltage (Vdc), across a single row of antidots, as a function 
of the driving amplitude with a frequency of 100 kHz. As 
explained in the text, three different dynamical regimes can 
be identified (labeled I, II and III) as a function of the ac am- 
plitude. The zoom in panel (e) gives the voltage dependence 
in region I. 



To simplify the presentation, we concentrate on the ex- 
perimental data obtained at the temperature T =1.430 
K and in a particular magnetic field if, which generates 
one flux quantum per each antidot. The equilibrium vor- 
tex configurations are visualized by plotting the square 
modulus of the order parameter, [Fig. [ijb)], and its 
real part. Re [Fig. [ijc)]. Since Re ip should change 
sign two times along a closed path encircling a single 
flux quantum, figure [TJc) indicates the presence of one 
^o-vortex trapped by each antidot. The stability of this 
vortex configuration leads to a minimum in dissipation 
(i.e. resistance) at this field [M], as shown in Fig. [ijd). 

Influence of the internal degrees of freedom .- Figure 
[TJe) shows the measured rectified (dc) voltage Vdc (pro- 
portional to the average vortex velocity) as a function 
of the amplitude of a sinusoidal excitation at a low fre- 
quency / = 100 kHz. The positive dc voltage sign cor- 
responds to particles moving preferentially in the —y di- 
rection. This motion, opposite to the preferential direc- 
tion at if = ifi, is caused by the competition between 
the vortex- vortex repulsion and the asymmetric pinning 
landscape ^ . A closer look at the amplitude dependence 
of Vdc reveals three distinct regions of rectification. At 
low ac amplitudes jac^ region I, a linear increase of Vdc is 
present (a close-up in the inset makes region I more vis- 
ible). With increasing ac amplitude this develops into a 



much steeper, non-linear increase in region II. The max- 
imum rectification accounts for the onset of region HI, 
where a sharp decrease of Vddjac) takes place. A single- 
particle model, as used in Ref. ^E], cannot describe the 
change in the rectifying properties between region I and 
II. As we demonstrate below, the sudden increase in the 
measured dc voltage results from a transformation of the 
moving vortices and their environment. 



0.4 



0.3 



> 

3 0.2 



0.1 



0.0 



— I \ 1 

1 (I) 1 


; 1. 1 

! (11) . 


It 

\k 

i ^ 


■ 

J I 


M 



(a) 



0.0 0.1 0.2 0.3 0.4 0.5 

/ (kA/cm') 



(b) 



4 

„ 2 

I ° 

> -2 
-4 



A. " ■ 


■A T ■ 


■ 1^ t ■ 1 ■ 1 ■ 


(I) 

> 1 - ' 


: ^ >!•■: 

(II) 

• i 


' 1 ' t 



0.0 2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 
ft ft ft 



FIG. 2. (a) Calculated dc voltage, Vdc, as a function of 
the amplitude of the ac current density for 7 = 10 and 
/ = 400 kHz. (b) Voltage drop between points (1) and (2) 
shown in Fig. 1(b), as a function of time for different ac cur- 
rents: (I) jac = 0.33 kA/cm^ (II) jac = 0.37 kA/cm^ (HI) 
jac = 0.39 kA/cm^. The dashed curve in panel (b) shows the 
dynamics of the applied current density. For the same ac cur- 
rents, snapshots of are shown in the insets of (a), labeled 
correspondingly (I), (II), (HI). The snapshots were taken dur- 
ing the ac cycle at maximum applied positive and negative 
currents. The thin and thick arrows indicate the direction of 
motion for Abrikosov and kinematic vortices, respectively. 

In order to elucidate the importance of the internal 
degrees of freedom of the vortices and the supercon- 
ducting condensate for the dc voltage read-out, we per- 
formed simulations using the generalized TDGL equa- 
tions [13]. An important model parameter, infiuenc- 
ing the internal structure of a moving vortex, is 7 = 
^kBTcTin\/u{l — T /Tc) / (ttK)^ with Tin, the inelastic col- 
lision time for electron-phonon scattering, and ^ 5.79. 
At 7 ^ 1 the healing time of the condensate is r ~ 
l\'^\^GL^ where tgl is the Ginzburg-Landau relaxation 
time [11] . In our calculations we have checked the behav- 
ior of the system for a wide range of 7 values. All other 
model parameters coincide with the experimental values. 

The dependence of Vdc on the driving amplitude jac is 
shown in Fig. j^Ja) for / = 400 kHz and 7 = 10. The 
TDGL simulations reproduce the qualitative features of 
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the experimental data [Fig. [TJe)], showing the same three 
distinct regions of rectification. Between 0.14 kA/cm^ 
and 0.37 kA/cm^ an increase of the dc voltage is ob- 
served with a discontinuity at 0.36 kA/cm^. At this point 
the dependence of the dc voltage on jac changes from 
an approximately linear increase (region I, 0.14 kA/cm^ 
< jac < 0.36 kA/cm^) to a steeper non-linear increase 
(region II, 0.36 kA/cm^ < j^^ < 0.37 kA/cm^). 

In order to reveal the origin of these features in Vddjac) 
we investigated the time evolution of the voltage differ- 
ence across a single row of antidots [between points (1) 
and (2) in Fig. ^h)]. At 0.33 kA/cm^ (region I) a clearly 
asymmetric motion of individual vortices is observed. A 
multitude of peaks in the voltage [Fig. [2|b)I] during one 
ac cycle indicates the passage of a series of vortices. A 
snapshot of [Fig. [2|a), Inset I] at maximum applied 
positive and negative currents indicates the presence of a 
vortex with a localized normal core at these low current 
drives. Since the dynamics of the condensate is limited by 
a specific time scale r, a drastic change occurs if the con- 
densate does not have enough time to fully recover before 
the passage of a new vortex [11]. Under these circum- 
stances the vortex will be deformed at higher velocities 
and a phase slip line (PSL) or kinematic vortex appears. 
A snapshot of \ip\'^ in region II at jac = 0.37 kA/cm^ 
[Fig. [2ja), Inset II] indeed shows the depletion of the or- 
der parameter along the direction of the vortex motion. 
In region II the broken symmetry of the vortex motion re- 
sults in the generation of such a PSL for positive applied 
currents only. The transformation between the motion of 
vortices with a localized core to a PSL exactly coincides 
with the boundary between the two different regimes of 
rectification. The motion of kinematic vortices is much 
faster and results in a highly increased voltage during 
the positive half-period of the ac drive [Fig. |2|b)II] and, 
as such, it explains the drastic increase of Vdc in region 
II. Above 0.37 kA/cm^ (region III) a PSL is created dur- 
ing both positive and negative half-periods of the applied 
current [panel (b)III and Inset III in Fig. [2], thus giving 
rise to a strong reduction of the average voltage signal at 
higher ac drives. It is important to note that in the PSL- 
like regime the vortex core elongates and it is impossible 
to see a vortex as a single entity. In this regime, the 
representation of this system by a point particle model 
in a periodic potential breaks down. Instead, a proper 
treatment of the vortex within the TDGL formalism is 
needed to explain the observed ratchet effects. 

Frequency dependence.- Since the transition from 
Abrikosov vortices to PSL is related to r, a change in 
the rectification properties is expected if the period of 
the driving force becomes comparable to this character- 
istic time scale. Unambiguous evidence of this change 
comes from Fig. |3] where we plotted the maximum recti- 
fied voltage, V^^^^\ as a function of the applied ac fre- 
quency (open circles). Above 300 kHz a strong decrease 
of rectification is observed until 3Mhz, beyond which only 



a small decay of the rectification signal is found. 

As discussed in the previous section, at low frequency 
drives, the maximum rectification is determined by the 
unidirectional creation of a PSL. The onset of PSL cor- 
responds to the regime in which the condensate is in- 
capable of full recovery in between two consecutive vor- 
tex passages. However, if the period of the driving force 
becomes comparable to the lifetime of a PSL, the PSL 
regime persists for both current directions. In the exper- 
iment this occurs above 3 Mhz. As a result, the PSL in 
this regime does not significantly contribute to the recti- 
fication signal (region II disappears) and a strong reduc- 
tion of the maximum rectification is expected. For those 
frequencies, the maximum of Ydc is achieved just before 
the onset of PSL, when rectification is dominated by the 
"conventional vortex hopping". A qualitatively similar 
frequency dependence is obtained in TDGL simulations 
[Inset of Fig. [s], indicating that the model captures the 
main physical mechanism. A comparison between the dif- 
ferent 7 values [7 = 10 (full circles) and 5 (full squares)] 
shows a shift of the overall behavior towards lower fre- 
quencies at higher 7 values. This effect is explained by 
the healing time of the condensate r being proportional 
to 7, so that the condensate acts slower for larger 7 val- 
ues. In line with this, the relatively low frequency scale 
in the experiment can be attributed to the fact that the 
7 value experimentally obtained for Al [15] is ~ 500. 
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FIG. 3. Maximum rectification (open circles) and the number 
of steps (open squares) as a function of the frequency / at 
T — 1.430 K and H — 2Hi. Inset: Maximum rectification, 
^i™^^^ (/iV), calculated within the TDGL model for 7 = 5 
(full squares) and 10 (full circles), as a function of / (MHz). 

Discrete rectification.- We concluded above that at suf- 
ficiently high frequencies the rectification is governed by 
the directional motion of Abrikosov vortices. In this 
regime a shorter timescale appears, associated with the 
traveling time of a vortex from cell to cell. The synchro- 
nization between an external excitation and the internal 
natural frequency should lead to a discrete hopping of 
vortices known as a vortex stepper motor. The aver- 
age voltage generated by a single moving vortex row is 
given by Vdc = ristep^of with ^0 = 2.07 x 10-^^ Tm^ 
the flux quantum, / the frequency in Hz, Ustep the num- 
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ber of steps taken by the vortex row. Knowing /, we 
can calculate the number of the unit cells the vortex is 
moving during a single ac cycle (as shown in Fig. [sj red 
circles). At low driving frequencies, this number is large 
{'^step > 100). As a result, any discreteness of the vortex 
motion becomes difficult to resolve. However, above 3 
MHz the number of steps taken per cycle is less than 10 
and discrete behavior of the vortex lattice should appear. 

Figure [4] shows a representative Vdc{Iac) curve (crosses) 
in the high frequency regime (/ = 6.5 MHz). Here a 
smoothed staircase structure is observed. The fact that 
the discrete hopping of vortices does not manifest itself as 
very sharp voltage steps, can be explained by the large 
number of antidot rows (333) which contribute to the 
measured signal. This could lead to statistical spread- 
ing. A common practice to unveil these blurred steps is 
to measure the derivative as shown in Fig. j4] (cir- 

cles). The minima in are used to determine the 

locations of the steps on the current axis and the height 
of the steps is defined as the voltage output exactly at 
this location. In the inset of Fig. [4] we plotted the height 
of first (squares) and second (triangles) voltage step to- 
gether with the theoretical prediction Vdc = nstep^of 
(red and black lines). The nice overlap between experi- 
ment and theory demonstrates that we indeed achieved a 
coherent stepper motor behavior of the whole Abrikosov 
lattice in the periodic pinning landscape. 
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FIG. 4. Vdc (crosses) and -^f^ (circles) as a function of lac 
with / = 6.5 MHz. Inset: Height of the first (squares) and 
second (triangles) voltage step as a function of /. The full 
lines indicate the expected values Vdc = ristep^of- 

To conclude, the main points of this letter are (i) the 
realization of the vortex stepper motor, particularly evi- 
dent at high frequency excitations, and (ii) the presence 
of a driving force window where the nature of the moving 
particle changes from an Abrikosov vortex to a phase slip 



line, when switching the direction of the driving force. 
The former effect is very sensitive to the particular value 
of the magnetic field, whereas the latter is a more robust 
effect. The formation of phase slip lines gives rise to the 
experimentally observed irreversibilities which have been 
previously explained by attributing an artificial mass to 
non-deformable vortices [7 . We believe that our results 
have relevance for other ratchet systems of deformable 
objects or media [4 . 
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